Role of microscopic phase separation in gelation of aqueous gelatin solutions 
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Using a unique home-made cell for four-contact impedance spectroscopy of conductive liquid 
samples, we establish the existence of two low frequency conductivity relaxations in aqueous solutions 
of gelatin, in both liquid and gel state. A comparison with diffusion measurements using pulsed field 
gradient NMR shows that the faster relaxation process is due to gelatin macromolecule self-diffusion. 
This single molecule diffusion is mostly insensitive to the macroscopic state of the sample, implying 
that the gelation of gelatin is not a percolative phenomenon, but is caused by aggregation of triple 
helices into a system-spanning fibre network. 



Gels are solutions where a small fraction of a gelling 
solute controls the phase of the whole system |l(. These 
gelling agents can produce macroscopic disordered net- 
works of interconnected molecules that span the whole 
solution. Physical properties of gels are mainly governed 
by the binding energy - usually of the order of the ther- 
mal energy - and the transition from a fluid (sol) state 
to a solid-like (gel) state is induced by lowering the tem- 
perature. The gelation temperature T g for this transition 
can be defined as the temperature where the sample stops 
to flow (typicaly around 25° C). This is correlated with 
the formation of the system-spanning network. 

In this Letter we study a low concentration aque- 
ous solutions of gelatin - a water soluble polypeptide^ 
made by denaturation from collagen, an irreversible pro- 
cess that destroys one collagen molecule to form three 
gelatin strands. At low temperatures (T < T g ) it forms 
macroscopically homogeneous and transparent thermore- 
versible (weak) gels (gel state) |3[. Above the gelation 
temperature (at T g ~ 40° C), the solution behaves as a 
normal viscous liquid (sol state). At the molecular level 
the polypeptide chains have a random coil conformation. 
By lowering the temperature they slowly change to a sin- 
gle helix conformation. If the concentration of gelatin 
is sufficient for overlap between nearby molecules, the 
formation of a triple-helix can be initiated by bonding 
with up to two other gelatin molecules. The resulting 
triple-helix is an insoluble semi-rigid rod with the par- 
tial structure of collagen due to incomplete renaturation 
but different energetics 3-6] . This process is the basis 
for two models of the system-spanning network genesis: 
(a) bond percolation[l[ 0) where macroscopic gel forms 
when the number of intermolecular bonds (triple helices) 
exceeds a treshold value and a single connected cluster 
is created; and (b) aggregation [7H9| where triple helices 
first aggregate into fibres, which then form a macroscopic 
fibre network. 

Although a widely studied system, there is no con- 
sensus on the nature of the sol-gel transition in gelatin. 
This is unsurprising because gelation is closely related 
to several other microscopic transitions that span the 
same temperature ranges like the coil-helix transition 



and the still not completely understood triple-helix for- 
mation. It is known that semi-rigid triple-helices con- 
tribute to the elasticity of gels, although it is still un- 
clear if through formation of a percolation cluster or a 
fibre network [lOj. Recent studies have shown that there 
is no simple connection between their macroscopic and 
microscopic properties ll|, |12[ , making the description of 
the gelatin gel as a percolation cluster questionable. 

Furthermore, the existence of a system-spanning net- 
work may impose a restriction on thermal concentra- 
tion fluctuations, implying an ergodic-nonergodic phase 
transition (dynamic arrest) [l3l similar to transi- 
tions in certain glassy systems [7J, |l5j. In our opin- 
ion the nonergodic nature of the transition masks its 
true behaviour [15j, causing inconsistencies in the results 
of light scattering diffusivity studies in the gel[I3, fisl — 
22j | . We propose that combined measurements of fre- 



quency dependent conductivity, and diffusivity of gelatin 
molecules by pulsed-field-gradient nuclear magnetic reso- 
nance (PFG-NMR) can resolve those difficulties. In par- 
ticular, we argue that the system-spanning network is 
formed by aggregation into fibre networks, as proposed 
for some biopolymer gelsQ. 

The sample used in this work is pig skin gelatin (Sigma- 
Aldrich, type A, Bloom 300). The solutions were pre- 
pared by dissolving gelatin in deionized water and stir- 
ring at 60°C for 2 hours, erasing any memory the system 
might have had. Sample volume was typically 75 [A, with 
gelatin concentrations c in the dilute gelling regimes be- 
tween the minimum gelling concentration c go i ~ 5 g/1 
and the overlap concentration c* ~ 25 g/1; and concen- 
trations corresponding to semi- dilute unentangled solu- 
tions (c* < c < c cn t ~ 150 g/1) (2^1 . Gelation is a slow 
kinetic process, meaning that it is impossible to stop the 
system at an intermediate step and measure the dynam- 
ics. Thus, we did all the measurements after a quench 
and wait time (non-equilibrium measurements[18]), such 
that we minimize the long length scale structural relax- 
ations but still obtain systematic gelation Tempera- 
ture measurements were done while adhering to a strict 
protocol: a quench (5 K/s) from 50°C to 5°C. After wait- 
ing for 1 hour for the gelatin to renaturate and relax0,0], 
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Figure 1. Conductivity spectra for a 85 g/1 gelatin solution 
at 42° C (circles) and 10° C (triangles). The real (filled sym- 
bols) and imaginary (empty symbols) part of the complex 
conductivity are shown. Solid lines are fits to two Cole-Cole 
relaxations [241]. Dotted lines show the individual contribu- 
tions of the two relaxation processes in Im(<r). 



Figure 2. Concentration scaling of the inverse characteristic 
time for the fast relaxation mode r/. The results are given 
for both the sol and gel state. Concentration c* = 25 g/1 
corresponds to the crossover between dilute and semi-dilute 
regimes. Inverse characteristic time scales as 1/t/ ~ c ' 3 in 
the gel sate and semi-dilute state, while in the dilute sol state 
it scales as 1/tj ~ c . 



the temperature was raised every 15 min and after equi- 
libration a spectrum was acquired. All temperature de- 
pendent data are, except when noted differently, given 
for 10°C (gel state) and for 42°C (sol state), well away 
from the temperature range needed for full denaturation 
or renaturation of the triple helices [H, 0). 

We measure impedance spectra of gelatin solutions in 
a planar four-contact arrangement, effectively eliminat- 
ing electrode polarization artefacts [2 5.]. This technique 
is similar to dielectric spectroscopy, with the difference 
that we detect the response of free charges, as opposed 
to bound charges (e.g. dipoles). Thus in our case the ap- 
propriate response function is conductivity a (uS) instead 
of permittivity e(uS). Fig. [1] shows (typical) measured 
spectra for 85 g/1 gelatin solution, for T — 10°C and 
42° C. We note two contributions to the conductivity re- 
laxation, the slow one and fast one. This is in contrast 
to other impedance spectroscopy results made with two- 
contact impedance cell (in frequency range from 1 Hz 
to 10 MHz) [2^, [27]], where low frequency contribution of 
the impedance is screened by the (unwanted) large elec- 
trode polarization effects, typical for conductive liquid 
samples. Other attempts to mitigate electrode polariza- 
tion effects had limited success (see references within Ref. 
[25[), proving that our custom made four-contact setup is 
an adequate tool for conductivity studies of liquid sam- 
ples at low frequencies. 

Gelatin can be considered a weak polyelectrolyte with 
a net negative charge of several percent of the number of 
monomeric units [28|. Our samples thus consist of free, 
weakly charged gelatin molecules, dissociated countcri- 
ons, and any excess salt contained in the sample. Mi- 
croscopically we can treat the large gelatin molecules 



as Brownian particles and use the corresponding equa- 
tion of motion (Langevin equation) to obtain a Drude- 
like relaxation of the frequency-dependent conductivity, 
a(uS) — Act/(1 + mjt)[2{J. This relaxation has the same 
form as a Debye relaxation in the dielectric constant. To 
include the possibility of a spread of characteristic time 
scales in each process, we analyse our data with the gen- 
eralized form of the Debye (Drude) formula, the so called 
Cole-Cole 24] functional form for complex frequency de- 
pendent conductivity, with fast and slow relaxation pro- 
cesses: 

-/ \ A07 Aa s . . 
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where r* s are the characteristic times, (3f. s exponents 
corresponding to the spread of relaxation times in the 
spectrum, Acr/ jS contributions of the (f)ast/(s)low pro- 
cess to the total conductivity, and a x the "infinite" fre- 
quency conductivity due to counterions. Full lines on 
Fig. [T] show Cole-Cole fits to real and imaginary parts 
of the measured conductivity at fixed temperature, and 
dotted lines show two individual contributions, 'slow'[3Cj 
and 'fast', to the Im(cr). From the fits, we can extract 
characteristic relaxation times t s ~ 1 s and 77 ~ 10 ms. 

In this work we focus on the faster of the two relax- 
ations (17 ~ 10 ms). Fig. [2] shows dependence of the 
fast process relaxation time, from Cole-Cole fits, on the 
gelatin concentration c. In dilute sol regime (c < c* 
and T = 42°C) relaxation time scales linearly with c _1 , 
while in semi-dilute sol regime (c* < c < c cn t) and all 
gel regimes (T = 10°C) we obtain an unusual scaling 
of Tf ~ c~ 03 . We also note that the relaxation is con- 
sistent with a Debye type relaxation, i.e. the exponent 
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f3f is always equal to 1 within the experimental error, 
implying no local heterogeneities for the charge carrier 
environment. The nature of the relaxation and the ob- 
tained exponents indicate that the fast relaxation can be 
interpreted as a conductivity contribution from gelatin 
molecule self-diffusion. 

To obtain additional information about this self- 
diffusion process we used pulsed field gradient (PFG) 
NM RI31II . with a standard stimulated echo sequence 
32L [33j. The technique enables direct measurement of 
the dynamic structure factor S (k,t) by spatially 'tag- 
ging' spins using magnetic field gradients (with a lin- 
ear relationship between wavenumber k and the gradi- 
ent strength), and recording their position after a time 
t. Experiments were performed on hydrogen nuclei at 
a NMR frequency of 360 MHz, in a home-made probe 
with built-in gradient coils. Dynamic structure factors 
S (k, t) obtained from PFG-NMR measurements contain 
two contributions: a fast 'normal diffusion' part at times 
~ f — 100 ms over long distances (~ 1 — 10 /an), con- 
sistent with the conventional time and k dependence 
S (k, t) ~ exp (— k 2 Dft/2>) and also observed in dynamic 
light scattering (DLS) studies [l7l [Lsl l3~ij] . and an anoma- 
lous, slow part at long times (~ 1 s). By comparison 
of the time-scales involved, we conclude that the two 
contributions correspond to the two relaxation processes 
observed in conductivity experiments. The effective dif- 
fusion coefficient for the fast process, Df, is obtained 
directly from the short-time slope in a In S vs. t plot. 
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Figure 3. a) The quantity nfqj obtained from diffusion 
and conductivity (see text for details). Lines are fits to 

n S q) ~ C °-98±0.03 for the SQl; and nfg 2 „ c 0.74±0.02 for ^ gd 

state, b) Temperature dependence of the normalized number 
of gelatin molecules lost to the gel macrostructure (rifq^ data 
for different concentrations were normalized to the same low 
and high temperature value). 



By combining the conductivity data and PFG-NMR 
diffusion we can obtain a relation between the concen- 
tration of charge carriers and dissolved gelatin. The 



Langevin equation approach provides us with a sim- 
ple (Einstein) equation which connects the character- 
istic relaxation height Act/ and the diffusion constant 
Df with the effective number of charge carriers, nf. 
A07 = nfqjDf/kT^^ (with q f the effective charge of 
the carriers). We can now use the concentration depen- 
dent values of Act/ and Df to obtain the concentration 
scaling of n/qi. This is shown in Fig. [3^ where a scaling 
of the form nfqj ~ c b is obtained with b = (0.98±0.03) in 
the sol and b = (0.74±0.02) in the gel state. This scaling 
is unchanged at the transition concentration c* between 
the dilute and semi-dilute regime, proving that we are in- 
deed dealing with self-diffusion of gelatin molecules. The 
sol state scaling shows that there is a linear relationship 
between the number of charge carriers and concentration 
of gelatin, while in the gel state we lose a concentra- 
tion dependent fraction of charge carriers due to gelation 
(see Fig. (3^,)- The reason for the incomplete renatura- 
tion most likely lies in the inherent randomness between 
gelatin molecules at the molecular level |3|. Due to a few 
orders of magnitude faster renaturation[35j and the lim- 
ited contribution of helix formation to rheologyjij we find 
that the n/q 2 ~ c~ a74 scaling in the gel state is consis- 
tent with the picture that a fraction of the molecules at 
any time are a part of the gel macrostructure and thus 
cannot contribute to the conductivity. Fig. [3Jd shows the 
reduced temperature dependence of n/cn normalized to 
their high temperature value, for several concentrations. 
All data conform to a master curve, with the point of 
inflection roughly corresponding to the gelation temper- 
ature. This is consistent with light polarization studies 
where the same behaviour is seen on the level of triple- 
helix formation (triple- helix melting curve) [J, |8| , |23| . l3C3j j . 
Preliminary circular dichroism (CD) results [30lj also in- 
dicate that a master curve exists for the concentration 
of triple-helices and that it has the same form, although 
the transition temperature obtained from conductivity is 
somewhat higher and with a different concentration de- 
pendence. The numerical values for nfqj for the case of 
~ 100 monomeric units in our gelatin indicate that about 
2-3% of the gelatin monomers have dissociated, as would 
be expected for a weak polyelectrolyte. 

We can also obtain a characteristic length scale L asso- 
ciated with the self-diffusion process by using the known 
scaling relation L 2 ~ DfTf, as shown in Fig. 2] The 
weak temperature dependence of L for two concentra- 
tions shown in the inset of Fig. H] shows that L is a well 
defined length scale, roughly independent of temperature 
and thus of the macroscopic state in the entire concen- 
tration range. However, the dynamics leading to a well- 
defined L is markedly different for dilute and semi-dilute 
solutions. The behaviour of dilute solutions in the sol 
state (c g < c < c* , T = 42°C) is consistent with a pic- 
ture of diffusing independent swollen coils: the scaling 
Tf ~ c _1 (Fig. [2]) shows that the time between collisions 
is proportional to the available volume per molecule, 
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Figure 4. Concentration scaling of the characteristic length 
scale for self-diffusion L 2 = DfTf. For concentrations larger 
than the overlap concentration c* = 25 g/l the data is con- 
sistent with a scaling of the form L 2 ~ c -1 ' 5 for both sol and 
gel. Below c* the exponent reduces to L 2 ~ c -0,85 for the 
sol state, while for the gel state scaling exponent is not well 
defined. Inset shows the characteristic length scale L as a 
function of temperature T. 



and for the characteristic length we obtain L 2 ~ c -0 85 , 
close to the scaling expected for non-overlapping spheres 
L 2 ~ c -2 / 3 . The value of the diffusion constant indicates 
a gyration radius for the swollen coils of R g « 20 nm. 
Scaling for the dilute gel state (c g < c < c* , T = 10°C) 
turns out not to be well defined, which we speculate is due 
to the complexity of the gel forming below the polymer 
overlap concentration, with possible competition between 
aggregation and percolation. 

In the semi-dilute sol and semi-dilute gel states (c > 
c*) we believe we are detecting translational diffusion re- 
sulting from reptation of gelatin macromolecules. We 
find L 2 ~ c~ 15 , consistent with the scaling of the mesh 
size from the De Gennes theory of polymers in good 
solvents [H, with the mesh size ~ 100 — 300 nm (slightly 
larger than indicated 

previously[I3,[Il|32i). DLS stud- 
ies report an analogous mode having a scaling in semi- 
dilute D ~ c _1 - 75 [l(| [13], while our results indicate a 
lower value D ~ c -12 . Although in the case of reptation 
one would expect a non-Debye relaxation and anoma- 
lous diffusion in PFG-NMR, we emphasize that we are 
working at longer length-scales than DLS: at our length 
scales (~ /im) the reptation dynamics renormalizes into 
the translational diffusion Df of the center of mass of 
the gelatin molecule [H| . This makes PFG NMR better 
suited for comparison with conductivity studies. Both 
the theory and our experimental data show that the cor- 
responding characteristic length L scales as the mesh size 
for (non-charged) semi-dilute polymers. If we now con- 
sider the gel state in the semi-dilute regime, we have the 
same length scaling and time dependence as in the sol. 
We thus conclude that the self-diffusion must be indepen- 



dent of the macroscopic state of the sample, depending 
only on the existence of a polypeptide mesh, originat- 
ing from either chain overlap or gelation. This confirms 
that rheological properties and diffusion of gelatin are 
independent; there are no heterogeneities in the local en- 
vironment of mobile gelatin molecules. We also observe 
that the length scale concentration dependence defining 
a mesh is almost unchanged due the sol-gel transition. 

The concentration of non-bound gelatin as implied by 
Fig. and 2] is sufficient to initiate the formation of 
triple-helices, but no such effect is seen in self-diffusion. 
This shows that after gelation we have two phases with 
a fast exchange of gelatin molecules, one analogous to 
the sol phase but less concentrated, and an aggregated 
phase of renaturated gelatin forming a fibre network. 
Gelatin, like collagen and similar systems [3^, can 
form fibres Q due to the attractive interaction between 
triple-helices in collagen [4l|. These interactions lead 
to spinodal decomposition in a bad solvent (e.g. water- 
methanol), but in pure aqueous dilute and semi-dilute 
solutions, like the one studied here, there is no compe- 
tition between gelation and precipitation because of the 



quality of solvent [42J, |43|. Thus the aggregation we ob- 
serve can be regarded as a form of microphase separation 
distinct from spinodal decomposition [![. The network of 
aggregated fibres is responsible for macroscopic proper- 
ties like elasticity, explaining why no simple correlation 
exists between gel strength (macrostructure) and diffu- 
sivity of probe particles [12j| and why microscopic probes 
like light scattering predict a different gelation transition 
point from macroscopic ones like rheology ll]. 

Because the gelation process was previously connected 
with percolation, current theories for similar systems do 
not account for the phase separation in the good solvent 



regime [44j,|45[ due to neglected microscopic cooperativity. 



Available theories for thermoreversible gelation (see ref- 
erences in 46 1 and 47-4t|) and numerical studies (Hoj do 
not predict the concentration scaling for the diffusion and 
relaxation times nor the gelling behaviour in the dilute 
regime. Our results lead to the conclusion that gelatin 
is closer to systems showing phase separation with gela- 
tion like colloid gelsfHlj and similar systems [H2I All 
this would make gelatin an excellent system for studying 
the complex phenomena of fibre networks 4^, 54],[55j|, e.g. 
glassy behaviour under stress [§]. 

In summary, by applying a four-electrode impedance 
spectroscopy method to dilute and semi-dilute gelling 
gelatin solutions we show the existence of two processes 
that contribute to the frequency dependent conductiv- 
ity. We identify the fast process as gelatin molecule 
self-diffusion due to long length scale reptation dynamics 
(translational diffusion). Self-diffusion at the crossover 
from dilute to semi-dilute state indicates that the mi- 
croscopic diffusivity of particles is decoupled from the 
macroscopic properties of the gel. Due to a loss of charge 
carriers in the gel we argue that upon gelation there ex- 
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ist two microphases in the system: the sol microphase 
where there is no interconnectivity and no difference from 
the sol macrophase, and the fibre microphase formed via 
aggregation of triple-helices into fibre networks. Our re- 
sults place semi-dilute (and dilute) gelatin solutions away 
from the standard bond percolation model, and they sug- 
gest that the gelation process goes via aggregation into 
fibre networks. This, although implied by the similar- 
ity of rheological data between known fibre-formers (e.g. 
actin)[56j], has not yet been established. Thus, studies 
of gelatin networks could give additional insight of the 
entire class of fibre forming systems. 
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